Uncertainty and Intelligence in Computational 
Stochastic Mechanics 


Bilal M. Ayyub 

Department of Civil Engineering 
University of Maryland 
College Park, MD 


127 




UNCERTAINTY AND INTELLIGENCE IN COMPUTATIONAL 
STOCHASTIC MECHANICS 


Bilal M. Ayyub 

Department of Civil Engineering 
University of Maryland 
College Park, MD 20742 


INTRODUCTION 

Classical structural reliability assessment techniques are based on precise and crisp (sharp) 
definitions of failure and non-failure (survival) of a structure in meeting a set of strength, function and 
serviceability criteria. These definitions are provided in the form of performance functions and limit 
state equations. Thus, the criteria provide a dichotomous definition of what real physical situations 
represent, in the form of abrupt change from structural survival to failure. However, based on 
observing the failure and survival of real structures according to the serviceability and strength criteria, 
the transition from a survival state to a failure state and from serviceability criteria to strength criteria are 
continuous and gradual rather than crisp and abrupt. That is, an entire spectrum of damage or failure 
levels (grades) is observed during the transition to total collapse. In the process, serviceability criteria 
are gradually violated with monotonically increasing level of violation, and progressively lead into the 
strength criteria violation. Classical structural reliability methods correctly and adequately include the 
ambiguity sources of uncertainty (physical randomness, statistical and modeling uncertainty) by 
varying amounts. However, they are unable to adequately incorporate the presence of a damage 
spectrum, and do not consider in their mathematical framework any sources of uncertainty of the 
vagueness type. Vagueness can be attributed to sources of fuzziness, uncleamess, indistinctiveness, 
sharplessness and grayness; whereas ambiguity can be attributed to nonspecificity, one-to-many 
relations, variety, generality, diversity and divergence. Using the nomenclature of structural reliability, 
vagueness and ambiguity can be accounted for in the form of realistic delineation of structural damage 
based on subjective judgment of engineers. For situations that require decisions under uncertainty with 
cost/benefit objectives, the risk of failure should depend on the underlying level of damage and the 
uncertainties associated with its definition. A mathematical model for structural reliability assessment 
that includes both ambiguity and vagueness types of uncertainty was suggested to result in the 
likelihood of failure over a damage spectrum. The resulting structural reliability estimates properly 
represent the continuous transition from serviceability to strength limit states over the ultimate time 
exposure of the structure. In this section, a structural reliability assessment method based on a fuzzy 
definition of failure is suggested to meet these practical needs. A failure definition can be developed to 
indicate the relationship between failure level and structural response. In this fuzzy model, a subjective 
index is introduced to represent all levels of damage (or failure). This index can be interpreted as either 
a measure of failure level or a measure of a degree of belief in the occurrence of some performance 
condition (e.g., failure). The index allows expressing the transition state between complete survival 
and complete failure for some structural response based on subjective evaluation and judgment. 
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STRUCTURAL RELIABILITY ASSESSMENT 


The reliability of an engineering system can be defined as its ability to fulfill its design purpose 
for some time period. The theory of probability provides the fundamental basis to measure this ability. 
The reliability of a structure can be viewed as the probability of its satisfactory performance according 
to some performance functions under specific service and extreme conditions within a stated time 
period. In estimating this probability, system uncertainties are modeled using random variables with 
mean values, variances, and probability distribution functions. Many methods have been proposed for 
structural reliability assessment purposes, such as First-Order Second Moment (FOSM) method. 
Advanced Second Moment (ASM) method, and computer simulation (Refs. 2 and 4). In this section, 
two probabilistic methods for reliability assessment are described. They are 1) advanced second 
moment (ASM) method, and 2) Monte Carlo Simulation (MCS) method with Variance Reduction 
Techniques (VRT) using Conditional Expectation (CE) and Antithetic Variates (AV). 

Advanced Second Moment (ASM) Method 

The reliability of a structure can be determined based on a performance function that can be 
expressed in terms of basic random variables X-s for relevant loads and structural strength. 
Mathematically, the performance function Z can be described as 

Z = Z(X„ X 2 ,..., X„J = Structural strength - load effect ( 1 ) 

where Z is called the performance function of interest. The failure surface (or the limit state ) of interest 
can be defined as Z = 0. Accordingly, when Z < 0, the structure is in the failure state, and when Z > 0 

it is in the safe state. If the joint probability density function for the basic random variables Xj's is 
/ = Z^ Xi Xn (x,, x 2 ,..., x„j , then the failure probability P f of a structure can be given by the integral 



/ m ,, 2 (*i, Xi,..., *„) dx l dx 2 ... dx n 


( 2 ) 


where the integration is performed over the region in which Z < 0. In general, the joint probability 
density function is unknown, and the integral is a formidable task. For practical purposes, alternate 
methods of evaluating P f are necessary. 

Reliability Index (Safety Index) 

Instead of using direct integration as given by Eq. 2, the performance function Z in Eq. 1 can be 
expanded using a Taylor series about the mean value of X's and then truncated at the linear terms. 
Therefore, the first-order approximate mean and variance of Z can be shown, respectively, as 


and 



(3) 

(4) 


130 



where Cov^X,, X ; j is the covariance of X, and X ; ; Z = mean of Z; and o 2 = variance of Z The partial 

derivatives of 3Z/3X, are evaluated at the mean values of the basic random variables. For statistically 
independent random variables, the variance expression can be simplified as 


CT 


2 

z 




(5) 


A measure of reliability can be estimated by introducing the reliability index or safety index P that is 
based on the mean and standard deviation of Z as 



(6) 


If Z is assumed to be normally distributed, then it can be shown that the failure probability P f is 


/» /= 1 -O(p) 


(7) 


where O = cumulative distribution function of standard normal variate. 

The aforementioned procedure of Eqs. 3 to 7 produces accurate results when the random 
variables are normally distributed and the performance function Z is linear. 

Nonlinear Performance Functions 

For nonlinear performance functions, the Taylor series expansion of Z is linearized at some 
point on the failure surface called design point or checking point or the most likely failure point rather 

than at the mean. Assuming the original basic variables X,'s are uncorrelated, the following 
transformation can be used: 


If X/s are correlated, they need to be transformed to uncorrelated random variables, as described by 

Thrift-Christensen and Baker (Ref. 33) or Ang and Tang (Ref. 2). The safety index P is defined as the 
shortest distance to the failure surface from the origin in the reduced T-coordinate system. The point 
on the failure surface that corresponds to the shortest distance is the most likely failure point. Using the 

original X-coordinate system, the safety index P and design point (xj, Xj, ..., X*j can be determined by 

solving the following system of nonlinear equations iteratively for P: 



f dz) 

W 

a *, 

n 

z 

/ = 1 

fdz ' 

w, 

2 

K- 

1/2 


X* = X, - a, Pa, 


(9) 


( 10 ) 
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z[x\,x\, ...,x;) = 0 


( 11 ) 


where a, = directional cosine; and the partial directives are evaluated at design point. Then, Eq. 7 can 
be used to evaluate P f . However, the above formulation is limited to normally distributed random 
variables. 

Equivalent Normal Distributions 

If a random variable X is not normally distributed, then it needs to be transformed to an 
equivalent normally distributed random variable. The parameters of the equivalent normal distribution 

X " and , can be estimated by imposing two conditions (Refs. 27 and 28). The cumulative 

distribution functions and probability density functions of a non-normal random variable and its 
equivalent normal variable should be equal at the design point on the failure surface. The first condition 
can be expressed as 


The second condition is 


0 


l < ) 



X , -X 


l = /M 


(12a) 


(12b) 


where F- t = non-normal cumulative distribution function; f = non-normal probability density function; 
O = cumulative distribution function of standard normal variate; and<j> = probability density function 
of standard normal variate. The standard deviation and mean of equivalent normal distributions can be 
shown, respectively, to be 


and 



(13) 




(14) 


Having determined c £ and X* for each random variable, P can be solved using the same procedure of 
Eqs. 9 to 11. 

The advanced second moment method is capable of dealing with nonlinear performance 
functions and non-normal probability distributions. However, the accuracy of the solution and the 
convergence of the procedure depends on the nonlinearity of the performance function in the vicinity of 
design point and the origin. If there are several local minimum distances to the origin, the solution 
process may not converge onto the global minimum. The probability of failure is calculated from the 

safety index P using Eq. 7 which is based on normally distributed performance functions. Therefore, 


132 



the resulting failure probability P f based on the ASM is approximate except for linear performance 
functions because it does not account for any nonlinearity in the performance functions. 


SOURCES AND TYPES OF UNCERTAINTY 

The following two viewgraphs show example sources of uncertainty, and a classification of 
uncertainty types. 


OBJECTIVES 

The objectives of this presentation were to generalize structural reliability assessment methods 
to account for ambiguity and vagueness sources of uncertainty, and demonstrate the developed methods 
using ship structures. A viewgraph is provided with a statement of objectives. 

Models and methods for merging different uncertainty sources in structural reliability 
assessment were described. The methods were presented in a finite element analysis framework. 

Also, intelligence in reliability computations with applications to marine vessels were discussed. 
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OBJECTIVES 


• Develop methods for structural reliability 
assessment based on a generalized treatment of 
uncertainty. 

• Define failure events over a damage spectrum. 

• Provide the reliability of the structure over the 
damage spectrum. 
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METHODOLOGY 


The following figure shows a procedure for an automated failure classification that can be 
implemented in a simulation algorithm for reliability assessment for ship structures as an example. The 
failure classification is based on matching a deformation or stress field with a record within a 
knowledge base of response and failure classes. In cases of no match, a list of approximate matches is 
provided, with assessed applicability factors. The user is prompted for any changes to the approximate 
matches and their applicability factors. In the case of a poor match, the user has the option of activating 
the failure recognition algorithm shown in the next figure to establish a new record in the knowledge 
base. The adaptive or neural nature of this algorithm allows the updating of the knowledge base of 
responses and failure classes. The failure recognition and classification algorithm shown in the figure 
evaluates the impact of the computed deformation or stress field on several systems of a structure. The 
impact assessment includes evaluating the remaining strength, stability, repair criticality, propulsion 
and power systems, combat systems, and hydrodynamic performance. The input of experts in ship 
performance is needed to make these evaluations using either numeric or linguistic measures. Then, the 
assessed impacts need to be aggregated and combined to obtain an overall failure recognition and 
classification within the established failure classes. The result of this process is then used to update the 
knowledge base. 

The development of a methodology for the reliability assessment of continuum ship structural 
components or systems requires the consideration of the following three components: (1) loads, (2) 
structural strength, and (3) methods of reliability analysis. Also, the reliability analysis requires 
knowing the probabilistic characteristics of the operational-sea profile of a ship, failure modes, and 
failure definitions. A reliability assessment methodology can be developed in the form of the following 
modules: operational-sea profile and loads; nonlinear structural analysis; extreme analysis and 
stochastic load combination; failure modes, their load effects, load combinations, and structural 
strength; library of probability distributions; reliability assessment methods; uncertainty modeling and 
analysis; failure definitions; and system analysis. Each module can be independently investigated and 
developed, although some knowledge about the details of other modules is needed for the development 
of a module. These modules are described by Ayyub, Beach and Packard (Ref. 6). 

Prediction of structural failure modes of continuum ship structural components or systems 
requires the use of nonlinear structural analysis. Therefore, failure definitions need to be expressed 
using deformations rather than forces or stresses. Also, the recognition and proper classification of 
failures based on a structural response within the simulation process need to be performed based on 
deformations. The process of failure classification and recognition needs to be automated in order to 
facilitate its use in a simulation algorithm for structural reliability assessment. The first figure shows a 
procedure for an automated failure classification that can be implemented in a simulation algorithm for 
reliability assessment. The failure classification is based on matching a deformation or stress field with 
a record within a knowledge base of response and failure classes. In cases of no match, a list of 
approximate matches is provided with assessed applicability factors. The user can then be prompted 
for any changes to the approximate matches and their applicability factors. In the case of a poor match, 
the user can have the option of activating the failure recognition algorithm shown in the second figure to 
establish a new record in the knowledge base. The adaptive or neural nature of this algorithm allows 
the updating of the knowledge base of responses and failure classes. The failure recognition and 
classification algorithm shown in the figure evaluates the impact of the computed deformation or stress 
field on several systems of a ship. The impact assessment includes evaluating the remaining strength, 
stability, repair criticality, propulsion and power systems, combat systems, and hydrodynamic 
performance. The input of experts in ship performance is needed to make these evaluations using either 
numeric or linguistic measures. Then, the assessed impacts need to be aggregated and combined to 
obtain an overall failure recognition and classification within the established failure classes. The result 
of this process is then used to update the knowledge base. 
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A prototype computational methodology for reliability assessment of continuum structures 
using finite element analysis with instability failure modes is described in this report. Examples were 
used to illustrate and test the methodology. Geometric and material uncertainties were considered in the 
finite element model. A computer program was developed to implement this methodology by 
integrating uncertainty formulations to create a finite element input file, and to conduct the reliability 
assessment on a machine level. A commercial finite element package was used as a basis for the 
strength assessment in the presented procedure. A parametric study for a stiffened panel strength was 
also carried out. The finite element model was based on the eight-node doubly curved shell element, 
which can provide the nonlinear behavior prediction of the stiffened panel. The mesh was designed to 
ensure the convergence of eigenvalue estimates. Failure modes were predicted on the basis of elastic 
nonlinear analysis using the finite element model. 

Reliability assessment was performed using Monte Carlo simulation with variance reduction 
techniques that consisted of the conditional expectation method. According to Monte Carlo methods, 
the applied load was randomly generated, finite element analysis was used to predict the response of the 
structure under the generated loads in the form of a deformation field. A crude simulation procedure 
can be applied to compare the response with a specified failure definition, and failures can then be 
counted. By repeating the simulation procedure several times, the failure probability according to the 
specified failure definition is estimated as the failure fraction of simulation repetitions. Alternatively, 
conditional expectation was used to estimate the failure probability in each simulation cycle in this 
study; then the average failure probability and its statistical error were computed. 

The developed method is expected to have significant impact on the reliability assessment of 
structural components and systems; more specifically, the safety and reliability evaluation of continuum 
structures, the formulation of associated design criteria, evaluation of important variables that influence 
failures, the possibility of revising some codes of practice, reducing the number of required costly 
experiments in structural testing, and the safety evaluation of existing structures for the purpose of life 
extension. The impact of this study can extend beyond structural reliability into the generalized field of 
engineering mechanics. 
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STRUCTURAL RELIABILITY ASSESSMENT 


The general performance function of a 
structural component or system according to a 
specified performance criterion is expressed as 
follows: 


Z = strength - load effect 
Z = g(X!, X 2 , X n ) 

where Xj = basic random variable 
g(*) > 0: survival event 
g(*) = 0: limit state 
g(*) < 0: failure event 
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The probability of failure is determined by 
solving the following integral: 


Pf= J I - I f x (Xi,X 2 ,-,X„) dxidx 2 -dx n 


where fx is the joint probability density function of 
X = {Xi,X 2 ,** , ,X n ) and the integration is performed over 
the range where g(*) < 0 
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UNCERTAINTIES 


I. Ambiguity: 


II. Vagueness: 


( 1 ) 

( 2 ) 

( 3 ) 

( 1 ) 

( 2 ) 


Physical randomness 
Statistical uncertainty 
Model uncertainty 

Definition of parameters 
Inter-relationships among the 
parameters 
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CRISP FAILURE MODEL 


Only two basic, mutually exclusive events, 
complete survival and complete failure , are 
considered, i.e., 

U -> {0, 1} 

where U = the universe of all possible outcomes 

0 = failure level of the event complete 

survival 

1 = failure level of the event complete 

failure 
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Failure Level, a 


complete 1.0 
failure 


complete 0.0 
survival 

Structural Response, R 
(e.g., curvature, deflection, etc.) 


Rf = structural response at failure 
R < Rf (a=0) : complete survival 
R = Rf : limit state 

R > Rf (a=l) : complete failure 
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FUZZY (CONTINUOUS) FAILURE MODEL 


A subjective index, failure level a, is introduced to 
represent the intermediate levels of damage, i.e., 

U A = { a : a e [0, 1] } 


where U = the universe of all possible outcomes 
a = 0 : complete survival 
0 < a < 1 : partial failure 
a = 1 : complete failure 


a can be interpreted as the degree of belief of a 
failure condition. 
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Failure Level, a 



Structural Response, R 
(e.g., curvature, deflection, etc.) 


Ri = lower bound of structural response 
R u = upper bound of structural response 
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Degree of Belief of an Event, a 



Structural Response, R 
(e.g., curvature, deflection, etc.) 


Event Number 

Definition 

1 

complete survival 

2 

low serviceability failure 

3 

serviceability failure 

4 

high serviceability failure 

5 

partial collapse 

6 

complete collapse 
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If definitions of failure events are interpreted as "at 
least low serviceability failure, serviceability 
failure, •••, or complete collapse," the above figure 
is modified to 


Degree of Belief of an Event, a 



Structural Response, R 
(e.g., curvature, deflection, etc.) 
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STRUCTURAL PERFORMANCE CURVE 



-0.8 -0.4 0.0 0.4 0. 


Curvature <|> ( x 10' 5 ) 




CRISP FAILURE MODEL 

FOR STOCHASTIC M-<P RELATIONSHIP 


Curvature (<J)) 



Moment (M) 
Load (L) 


fL = probability density function (pdf) of L 
fjvi = conditional pdf of M at = <|)f 
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The probability of failure is evaluated as 

Pf = Prob { a = 1 } 

= Prob { L > (M at <|> = <t>f) } 

oo 

= J Prob{L > (m at([)=<|)f)} fM(m) d m 
o 

oo 

= J {l-F L (m))f M (m)dm 
0 

where Fl = cumulative distribution function of L 
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FUZZY FAILURE MODEL 

FOR STOCHASTIC M-0> RELATIONSHIP 


Curvature (())) 



fL = probability density function (pdf) of L 
fM = conditional pdf of M at = c()f 
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The probability of failure is evaluated as 

Pf (0Cf) = Prob { Zlcc=ccf < 0 } 

= Prob { M(<|) f ) - L < 0 } 

oo 

= J Prob{L > (m at(()=c()f)} fM(m) d m 
o 

oo 

= J {l-F L (m)}fM(m)din 
0 

where Fl = cumulative distribution function of L 
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AVERAGE PROBABILITY OF FAILURE 


I. Crisp Failure Model: 

Pf, avg = Pf 

II. Fuzzy Failure Model: 

• Arithmetic average: 

1 


J Pf(a) da 



J da 
o 


• Geometric average: 

1 

J logio ( p f(a)) da 
o 

logio(Pfg) = 1 

J da 
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EXAMPLE I 

Consider the following performance function: 
Z = M - L = M(<|0 - L 

where M = resisting moment (ft-tons) 

L = external load (ft-tons) 
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• Crisp Failure Model 

The curvature at failure is specified as 
<|>f = 0.30 x 10-5 

The statistical characteristics of external moment 
(L) and resisting moment (Mf): 


Random 

Variable 

Mean 

Value 

Coefficient 
of Variation 
(COV) 

Probability 
Distribution Type 

L 

100x103 

ft-tons 

0.30 

Extreme Value Type I 

M f 

244x103 

ft-tons 

0.10 

Normal 
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Fuzzy Failure Model 


Failure Level (a) 



0.2 0.25 0.275 0.325 0.35 0.4 

-5 

Curvature <|> ( x 10 ) 




Probability of Failure ( xlO * 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Failure Level ( a ) 


The values of Pf were calculated using 1000 simulation cycles. 
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Probability of Failure ( xlO ‘ 



The values of Pf were calculated using 1000 simulation cycles. 
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Probability of Failure 


0.005 


0.004 


0.003 

0.002 


0.001 


0.000 



0 200 400 600 800 1000 1200 


Number of Simulation Cycles 
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COV(P 
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0.25 
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Average Probability of Failure for Example I 


Failure 

Model 

Curvature at Failure 

4 > f 

Arithmetic 
Average of 
Probability 
of Failure 

Geometric 
Average of 
Probability 
of Failure 

Fuzzy 1 

0.2xl0* 5 to 0.4x10*5 

9.137x10*3 

2.320x10-3 

Fuzzy 2 

0.25x10*5 to 0.35x10*5 

2.752x10*3 

1.851x10*3 

Fuzzy 3 

0.275x10*5 to 0.325x10*5 

2.086x10*3 

1.854x10-3 

Crisp 

0.3x10*5 

1.973x10-3 

1.973x10*3 
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EXAMPLE II - CASE A 



160 


<j> ( X 10 *5) 




Degree of Belief ( a ) 



.0001 



Event 1 
Event 2 
Event 3 
Event 4 
Event 5 
Event 6 


Probability of Failure Occurrence 


The values of P f were calculated using 1000 simulation cycles. 
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Average Probability of Occurrence - Case A 


Event 

No. 

Definition 

Arithmetic Average 
of Probability 
of Occurrence 

Geometric Average 
of Probability 
of Occurrence 

1 

complete survival 

0.940 

0.940 

2 

low serviceability 
failure 

1.611x10-2 

1.583x10-2 

3 

serviceability 

failure 

8.718x10-3 

8.396x10-3 

4 

high serviceability 
failure 

4.944x10*4 

4.782x10-4 

5 

partial collapse 

1.439x10-4 

1.424x10*4 

6 

complete collapse 

2.847x10-4 

2.846x10*4 
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EXAMPLE II - CASE B 
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Degree of Belief ( a ) 




Event 1 
Event 2 
Event 3 
Event 4 
Event 5 
Event 6 


Probability of Failure Occurrence 


The values of Pf were calculated using 1000 simulation cycles. 
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Average Probability of Occurrence - Case B 


Event 

No. 

Definition 

Arithmetic Average 
of Probability of 
Occurrence 

Geometric Average 
of Probability of 
Occurrence 

1 

complete survival 

0.940 

0.940 

2 

low serviceability 
failure 

2.619x10-2 

2.555x10-2 

3 

serviceability 

failure 

1.008x10-2 

9.722x10-3 

4 

high serviceability 
failure 

9.388x10-4 

9.206x10-4 

5 

partial collapse 

4.444x10-4 

4.424x10-4 

6 

complete collapse 

2.847x10-4 

2.846x10-4 
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UNCERTAINTY MEASURES 

• Hartley Measure: Set theory 

• Shannon Entropy: Probability theory 

• Measure of Fuzziness: Fuzzy set theory 

• U-Uncertainty: Possibility theory 

• Measure of Dissonance: Theory of evidence 

• Measure of Confusion: Theory of evidence 
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UNCERTAINTY MEASURES 



A 

a 

c 

1 D 

E 

1 F 

G 

1 

Uncertainty measure 

Type of uncertainty 

Type of sets or events 

Hicsxnn 

Comments 


Reference 

2 

Hartley 

ambiguity 

crisp 

set 

A basic discrete measure. 

A larger number of outcomes 
means larger uncertainty 

[0.-) 

Hartley (1928) 

n 

4 

5 

6 

Shannon Entropy 

ambiguity 

crisp 

set 

and 

probability 

The closer the outcomes 
to an equal liklihood, the 
larger the uncertainty 

(0.-) 

Shannon (1946) 

■■ 

mm 

Dl 

□ -uncertainty 

ambiguity 

crisp 

set 

and 

possibility 

Possibiiistic counterpart 
to Shannon entropy and 
generalization to Hartley 
measure 

(0.-) 

Higashi and Klir [1903] 

BQ 

m 

KB 

m 

Fuzziness measure 

vagueness 

and 

ambiguity 

fuzzy 

set 

and 

fuzziness 

Measures the lack ot 
distinction between a set 
and its complement 

[0.-) 

Deluca and Termini 
[1972.1974.1 977] 

KD 

ID 


— 

crisp 

set 

and 

evidence 

Measures conflict ol 
evidence using theory of 
evidence 

10.-) 

Yager [1983] 

ID 

m 

ED 


BUSH 

crisp 

set 

and 

evidence 

Measures confusion of 
evidence using theory of 
evidence 

[0.-) 

Hohle [1981] 

ID 

ED 
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Failure Classification 
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Failure Recognition and Classification 
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Finite Element Mesh of the Stiffened Panel 
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Web Height Variability and Web Tilting 






VS EB BOW INC {LEVELS - 0 , I , ■ » d 2) 


TEN GENERATED RANDOM VARIABLES: X BOl, X BOl, X BOi, X BO<, 


X B 0 5 


Web Bowing 






Flange Width Variability and Range Tilting 
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Methodology for ith Simulation Cycle 











Write 

► strength and 
probability of failure 


Update files 
(SUMPf, SUMPfZ) 


C-Shell Script Flow Chart for ith Simulation Cycle-Machine Level 





















Geometric And Material Random Variables for the Stiffened Panel 


Variable 

no 

Geometrical variables 

Notation 

Mean value 

Coefficient of 
variation 
(COV) 

Standard 

deviation 

1 

Plate size (mm) 

LOi 

854 


4.0 

2 

Plate thickness (mm) 

to 

3.0 

4% 

0.12 

3 

Web thickness (mm) 

ti 

4.9 

4% 

0.196 

4 

Flange thickness (mm) 

t 2 

5.84 

4% 

0.234 

5 

Plate-out of plane distortion 
(mm) 

zPOi 

0.0 


1.0 

6 

Web height (mm) 

Llji 

31.08 

2.5% 

0.77 

7 

Web tilting (mm) 

X2ji 

0.0 


0.5 

8 

Web bowing (mm) 

XBOi 

0.0 


0.1 

9 

Flange width (mm) 

L2ji 

25.4 

2.5%' 

0.635 

10 

Flange tilting (mm) 

Z2i0, Z2iL 

0.0 


0.2 

11 

Modulus of elasticity (MPa) 

E 


4% 

8320 

mm 

Poisson's ratio 

V 




13 

i 

Yield stress (KPa) 1 

Fy 

250000 

7% 

17500 


1 Nominal yield stress = 240000 kPa 
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Thicknesses and Plate Geometric Variables 


Variable no. 

Geometrical variables 

Notation 

Mean value 
(mm) 

Coefficient 
of variation 
(COV) 

Standard 

deviation 

(mm) 

global 

local 

1 

1 

Panel width (side 1 ) 

Lo. 

854.0 


4.0 

2 

2 

Panel width (side 3) 

L-03 

854.0 


4.0 








3 

3 

Plate (hickness 

*P 

3.0 

4% 

0.12 

4 

4 

Web thickness 

tw 

4.9 

4% 

0.196 

5 

5 

Flange thickness 

tf 

5.84 

4% 

0.234 








6 

6 

Plate-out of plane distortion 
(corner2) 

Zpo2 

0.0 


1.0 

7 

7 

Plate-out of plane distortion 
(corner3) 

Zp<)3 

0.0 


1.0 

8 

8 

Plate-out of plane distortion 
(corner4) 

Zf04 

0.0 


1.0 


180 




Web Height Variables 


Variable no. 

Geometrical variables 

Notation 

Mean value 
(mm) 

Coefficient of 
variation 
(COV) 

Standard 
deviation (mm) 




i 

Height of web no. 1 (side 1) 

Li n 

31.08 

2.5% 

0.77 


2 

Height of web no. 2 (side 1) 

Lli2 

3 1 .08 

2.5% 

0.77 

B 

3 

Height of web no. 3 (side 1 ) 

Lj|3 

31.08 

2.5% 

0.77 

B 

B 

Height of web no 4 (side 1) 

Liu 

31.08 

2.5% 

0.77 

B 

B 

Height of web no. 5 (side 1 ) 

L||5 

31.08 

2.5% 

0.77 


Hi 






15 

B 

Height of web no 1 (side 3) 

■ 

— 

2.5% 

0.77 

B 

7 

Height of web no. 2 (side 3) 

L|32 

31.08 

2.5% 

0.77 

B 

8 

Height of web no. 3 (side 3) 

Lm 

31.08 

2.5% 

0.77 

18 

9 

Height of web no. 4 (side 3) 

L| 34 

31.08 

2.5% 

0.77 

19 

10 

Height of web no. 5 (side 3) 

Li 15 

31.08 

2.5% 

0.77 
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Web Tilting Variables 


Variable no. 


global local 


20 


Geometrical variables 


Tilting of web no.l (side 1) 


Tilting of web no. 2 (side 1 ) 


Tilting of web no. 3 (side 1) 


Tilting of web no.4 (side 1) 


Tilting of web no.5 (side 1) 


Notation Mean value Coefficient of Standard deviation 
(mm) variation (mm) 

(COV) 


Tilting of web no.l (side 3) 


Tilting of web no. 2 (side 3) 


Tilting of web no. 3 (side 3) 


Tilting of web no.4 (side 3) 


Tilting of web no.5 (side 3) 


X 23 I 

0.0 

X232 

0.0 

X233 

0.0 

X234 

0.0 

X235 

0.0 


Variable no. 


Web Bowing Variables 


Geometrical variables 


Notation 


Mean value 
(mm) 


Coefficient Standard 

of variation deviation (mm) 
(COV) 


Bowing of web no. I (side 3 ) X B oi 


Bowing of web no. 2 (side 3) X B o 2 


Bowing of web no. 3 (side 3) X B o3 


Bowing of web no. 4 (side 3) X B o4 


Bowing of web no. 5 (side 3) 



































































Flange Width Variables 


Variable no. 

Geometrical variables 

Notation 

Mean value 
(mm) 

Coefficient of 
variation 
(COV) 

Standard 
deviation (mm) 


local 



Hi 

Width of flange no. 1 (side 1) 

mm 

25.4 

2.5% 

0.635 

31 

2 

Width of flange no. 2 ( side 1) 

mm 

25.4 

2.5% 

0.635 

32 

3 

Width of flange no. 3 ( side 1) 

mm 

25.4 

2.5% 

0.635 

33 

4 

Width of flange no. 4 ( side 1) 


25.4 

2.5% 

0.635 

34 

5 

Width of flange no. 5 ( side 1) 

— 

25.4 

2.5% 

0.635 








35 

6 

Width of flange no. 1 ( side 3) 

f-23l 

25.4 

2.5% 

0.635 

36 

7 

Width of flange no. 2 ( side 3) 

1-232 

25.4 

2.5% 

0.635 


8 

Width of flange no. 3 ( side 3) 


25.4 

2.5% 

0.635 

MM 

9 

Width of flange no.4 ( side 3) 


25.4 

2.5% 

0.635 

39 

10 

Width of flange no.5 ( side 3) 

n 

25.4 

2.5% 

0.635 
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Flange Tilting Variables 


Variable no. 

Geometrical variables 

Notation 

Mean value 
(mm) 

Coefficient of 
variation 
(COV) 

Standard 

deviation 

(mm) 

global 

local 

40 

1 

Tilting of flange no.l (side 1) 

Z210 

0.0 


0.2 

41 

2 

Tilting of flange no. 2 (side 1) 

Z220 

0.0 


0.2 

42 

3 

Tilting of flange no. 3 (side 1) 

Z230 

0.0 


0.2 

43 

4 

Tilting of flange no. 4 (side 1) 

Zwo 

0.0 


0.2 

44 

5 

Tilting of flange no.5 (side 1) 

Z250 

0.0 


0.2 








45 

6 

Tilting of flange no.l (side 3) 

Z21L 

0.0 


0.2 

46 

7 

Tilting of flange no.2 (side 3) 

Z 22 L 

0.0 


0.2 

47 

8 

Tilting of flange no. 3 (side 3) 

Z 23 L 

0.0 


0.2 

48 

9 

Tilting of flange no.4 (side 3) 

Z 24 L 

0.0 


0.2 

49 

10 

Tilting of flange no.5 (side 3) 

Zm 

0.0 


0.2 
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Material Variability 


Variable no. 


global local 





Material variables 


Modulus of elasticity of 
plate material (MPa) 


Notation 


2 Modulus of elasticity of web 
material (MPa) 


52 3 Modulus of elasticity of 

flange material (MPa) 


Poisson s ratio of plate 


Poisson’s ratio of web 


55 6 Poisson's ratio of flange 


Yield stress of plate (kPa) 1 


57 8 Yield stress of web (kPa) 1 


| 9 Yield stress of flange (kPa) 1 

Nominal yield stress = 240000 kPa 



Coefficient of Standard 

variation (COV) deviation 




250000 

7% 

17500 

250000 

7% 

17500 

250000 

1 % 

17500 
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Buckling Shape of the Stiffened Panel 


186 





Axial Strength Stati: 



Normalized Axial Strength-500 cycles, Average = 1.002 



of the Stiffened Panel-500 Cycles 



































Stiffened Panel Subjected to Concentric Axial Compression and Lateral Pressure 
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PARAMETRIC ANALYSIS 


A parametric analysis was conducted for the axial strength and failure probability of the panel. 
The analysis was carried out by individually varying the coefficients of variation or standard deviations 
of the basic random variables. The notations, mean values, and ranges of COV and standard deviations 
of the random variables are given in the following table. The following observations were developed 
based on the results of the parametric analysis using 100 simulation cycles: 

• For the plate width, a figure shows that increasing the COV from 0.47% to 0.94%, the normalized 
strength decreases from 1.007 to 0.988, the COV of the axial strength decreases from 8.87% to 

7.79%, and the average of probability of failure decreases from 8.20 x 10' 4 to 6.45 x 10' 4 . 

• For the plate-out of plane distortion, a figure shows that increasing the standard deviation from 1 .0 to 
3.0, the normalized strength increases from 1.007 to 1.009, the COV of the axial strength decreases 

from 8.87% to 7.42%, and the average of probability of failure decreases from 8.20 x IO* 4 to 6.45 x 

10- 4 . 

• For the web height, a figure shows that increasing the COV from 2.5% to 5.0%, the normalized 
strength increases from 1.007 to 1.01 1, the COV of the axial strength decreases from 8.87% to 

7.37%, and the average of probability of failure decreases from 8.20 x IO -4 to 3.94 x IO -4 . 

• For the web tilting, a figure shows that increasing the standard deviation from 0.2 mm to 0.5 mm, the 
normalized strength increases from 1.005 to 1.007, the COV of the axial strength increases from 

7.17% to 9%, and the average of probability of failure increases from 5.0 x 10 -4 to 8.45 x 10‘ 4 . 

• For the web bowing, a figure shows that increasing the standard deviation from 0. 1 mm to 0.2 mm, 
the normalized strength decreases from 1 .007 to 0.99, the COV of the axial strength decreases from 

9.0% to 7.8%, and the average of probability of failure decreases from 8.20 x IO -4 to 4.84 x 10" 4 . 

• For the flange width, a figure shows that increasing the COV from 2.5% to 5.0%, the normalized 
strength decreases from 1.007 to 1.004, the COV of the axial strength decreases from 9.0% to 

7.26%, and the average of probability of failure decreases from 8.20 x IO -4 to 2.23 x 10' 4 . 

• For the flange tilting, a figure shows that increasing the standard deviation from 0.2 mm to 0.5 mm, 
the normalized strength decreases from 1 .007 to 0.995, the COV of the axial strength decreases from 

9.0% to 8.0%, and the average of failure probability increases from 8.20 x IO 4 to 1.77 x 10 -3 . 

• For the thicknesses, a figure shows that increasing the COV from 4.0% to 8.0%, the normalized 
strength decreases from 1.007 to 0.994, the COV of the axial strength increases from 8.87% to 

13.0%, and the average of probability of failure decreases from 8.28 x IO -4 to 1.53 x 10 2 . 

• For the modulus of elasticity, a figure shows that increasing the COV from 4.0% to 8.0%, the 
normalized strength decreases from 1 .007 to 1 .002, the COV of the axial strength remains constant at 

the value of 8.90%, and the average of probability of failure decreases from 8.20 x 10 -4 to 1.49 x 
IO 3 . 


The above failure probability observations were based on results from 100 simulation cycles. 
The number of simulation cycles might not be adequate for obtaining accurate failure probability 
results, but it is sufficient for determining the axial strength. The number of cycles was limited to 100 
in order to make the study feasible within the planned time frame of the project. 
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A table shows a summary of the results of the parametric study. According to the table, 
variations in the variability of plate size and web bowing produced the largest effect on the mean axial 
strength ratio; whereas variations in the variability of thicknesses of the plate, webs, and flanges 
produced the largest effect on the coefficient of variation of the axial strength. 
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Variation of Coefficient of Variation or Standard Deviation 


Variable 

no. 

Geometrical variables 

Notation 

Mean value 

Coefficient of 
variation 
(COV) 

Standard deviation 

1 

Plate size (mm) 

LOi 

854 


4.0 to 8 

2 

Plate thickness (mm) 

to 

3.0 

4 to 8% 


3 

Web thickness (mm) 

ti 

4.9 

4 to 8% 


4 

Flange thickness (mm) 

t 2 

5.84 

4 to 8% 


5 

Plate-out of plane 
distortion (mm) 

zPOi 

0.0 


1.0 to 3.0 

6 

Web height (mm) 

Llji 

31.08 

2.5 to 5% 


7 

Web tilting (mm) 

X2ji 

0.0 


0.2 to 0.5 

8 

Web bowing (mm) 

XBOi 

0.0 


0.1 to 0.2 

9 

Flange width (mm) 

L2ji 

25.4 

2.5 to 5% 


10 

Flange tilting (mm) 

Z2i0, Z2iL 

0.0 


0.2 to 0.5 

11 

Modulus of elasticity 
(MPa) 

E 

208000 

4 to 8% 
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Parametric Analysis Results 


Variable 

no. 

Geometrical Variables 

Mean 

value 

Variation 

of 

coefficient 
of variation 

Variation of 
standard 
deviation 

Effect on axial 
strength ratio 

Effect on 
coefficient of 
variation of 
strength 

1 

Plate size (mm) 

854 




Medium/Low 

2 

Plate thickness (mm) 

3.0 



Medium 

High 

3 

Web thickness (mm) 

mm 

4 to 8% 

0.196 to 0.392 

Medium 

High 

4 

Flange thickness (mm) 

5.84 

4 to 8% 

0.234 to 0.468 

Medium 

High 

5 

Plate-out of plane 
distortion (mm) 

0.0 


1.0 to 3.0 

Low 

Medium/Low 

6 

Web height (mm) 

31.08 

2.5 to 5% 

0.77 to 1.54 

Low 

Medium/Low 

7 

Web tilting (mm) 

0.0 


0.2 to 0.5 

Low 

Medium/Low 

8 

Web bowing (mm) 

0.0 


0.1 to 0.2 

High 

Medium/Low 

9 

Flange width (mm) 

25.4 

2.5 to 5% 

0.635 to 1.27 

Low 

Medium/Low 



0.0 


0.2 to 0.5 

Medium 

Medium/Low 

11 

Modulus of elasticity 
(MPa) 

208000 

4 to 8% 

8320 to 16640 

Low 

None 
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RECOMMENDATIONS FOR FUTURE WORK 


Based on this study, the following recommendations for future work are provided: 

The feasibility of using the developed method for complex structures with multiple failure modes 
needs to be investigated. The structures need to be selected such that methods for failure recognition 
and classification as previously demonstrated can be developed. 

The effects of failure recognition and classification for continuum structures on reliability estimates 
need to be studied. 
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